Majorana modes and complex band structure of quantum wires 
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We describe Majorana edge states of a semi-infinite wire using the complex band structure ap- 
proach. In this method the edge state at a given energy is built as a superposition of evanescent 
waves. It is shown that the superposition can not always satisfy the required boundary condition, 
thus restricting the existence of edge modes. We discuss purely ID and 2D systems, focussing in 
the latter case on the effect of the Rashba mixing term. 

PACS numbers: 73.63.NM,74.45.+c 



I. INTRODUCTION 

The realization of Majorana states as specific exci- 
tations of semiconductor quantum wires and of other 
condensed matter systems has attracted much interest 
recently*^— Experimental evidences of these peculiar 
states in quantum wires and topolo gica l insulators have 
been presented in Refs. [23l - l25j and [26. 27}, respectively. 
To assess these evidences, theoretical works are presently 
focussing on a deeper clarification of the physical prop- 
erties of Majorana states, including analysis of stabil- 
ity limitations and robustness against disorder and other 
perturbations^ — 

As most characteristic properties, a Majorana state 
is a zero-energy mode, thus degenerate with the system 
ground state, that is localized on the system edges or 
interfaces. Additionally, an energy gap with neighbor- 
ing states protects the Majorana modes against decoher- 
ence. In semiconductor wires, the subject of this work, 
such zero-energy modes can be induced by the combined 
action of the following three mechanisms: superconduc- 
tivity, Rashba spin-orbit coupling and Zeeman magnetic 
field. Superconductivity introduces particle-hole dual- 
ity, with excitations lying symmetrically at positive and 
negative energies with respect to the chemical potential. 
Rashba spin-orbit coupling introduces helicity by associ- 
ating motion in space with spin. Thirdly, the Zeeman 
coupling with an applied magnetic field can be viewed 
as a tunable interaction with which the system can be 
placed in different regimes. 

The edge character of the Majorana modes suggests a 
link with the physics of surfaces. Indeed, the similarities 
with the Schockley states of finite crystals^ have been 
pointed out in Ref. [22| • We explore in this work another 
link with surface physics, namely the description of Ma- 
jorana modes using the complex band structure of the 
system. The complex band structure of crystals (or peri- 
odic systems in general) is a well established theoretical 
formalism allowing the generalization from bulk propa- 
gating states to states that are localized around defects 
or interfaces of the crystal^ The underlying idea of the 
generalization is the following: while usual propagating 
states are described by sustained waves with a real wave 



number k, there are other solutions with complex wave 
numbers representing evanescent waves whose amplitude 
decreases in a given direction (increases in the reversed 
one). 

Waves with complex k cannot be physically realized 
for infinitely long distances from a given point in both 
directions due to the divergent behavior in one of them; 
but they can in a finite domain, or even in a semi-infinite 
one provided the infinite distance is in the direction of 
decaying amplitude. The latter situation is precisely the 
one we expect for edge Majorana modes. While for prop- 
agating states the analysis is usually done in terms of the 
band structure {E a (k)}, where a is the band index and 
k is an arbitrary real wave number, the complex band 
structure is presented as {k a (E)}, where for any given 
energy E there is a set of allowed complex wave numbers 
k a . Of course, if the energy does not allow propagat- 
ing states all k a will have an imaginary part. In general, 
however, propagating and evanescent modes may coexist. 

In this work we determine numerically the complex 
band structure of a semiconductor quantum wire with the 
above mentioned interactions. Previously^! we studied 
the complex band structure in absence of superconduc- 
tivity and Zeeman field. Our present numerical approach 
to Majorana physics complements the analytical meth- 
ods used in Ref. [2l[. We suggest robust numerical algo- 
rithms to obtain the complex band structure of ID and 
2D quantum wires in general, discussing the conditions 
for the existence of Majorana modes in the semi-infinite 
system. In 2D we focus on the role of the Rashba mix- 
ing, giving evidence that the Rashba mixing mechanism 
hinders the coexistence of two Majorana modes. This is 
in agreement with the findings of numerical diagonaliza- 
tions for the finite system.— Our results correspond to a 
continuum Hamiltonian of quantum wires and, from this 
point of view, they also complement approaches based on 
tight binding models^^ 

The analysis in terms of the complex band structure 
gives the precise conditions for the existence of zero- 
energy edge modes, as an alternative to the analysis 
based on propagating bands. In particular, it clarifies 
why the topological phase transition from absence to 
presence of a Majorana mode is signaled by the van- 
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ishing and reopening of the gap for propagating states 
with k — 0, as discussed in Ref. 0- Sections II and III 
are devoted to the ID and 2D cases, respectively. The 
first one is meant as a simple yet illustrative case of the 
more realistic situation including transverse degrees of 
freedom. 



II. THE ID CASE 

A. Hamiltonian and complex band structure 

Let us first consider a purely ID model as in Ref. 2l|, 
with motion constrained to be on the x axis. The Hamil- 
tonian reads in this case 



n 1D = 

1 2m 



A B (T • n + A t x 



a 

—p x a y T z 



(1) 



where the Pauli operators for spin and isospin (in 
electron-hole space) are represented by the cr's and r's, 
respectively. The successive energy contributions to Eq. 
(fTJ) are: kinetic, Zeeman (depending on parameter As), 
pairing (parameter Ao) and Rashba spin-orbit coupling 
(parameter a). In addition, Eq. (TTJ) includes the chemical 
potential (/x) and the magnetic field orientation (n). The 
system exhibits very different physics for parallel (n — x) 
and transverse (n = y, z) field orientations. Since in ID 
all transverse orientations are equivalent we will consider 
a field lying in the xy plane with an azimuthal angle 
< 4>b < 2tt. 

The band structure of the infinite system is obtained 
from Schrodinger equation 



(2) 



where the wave function variables are the space coordi- 
nate x £ (-co, +oo), the spin rj a 6 {t,i} and the isospin 
i] r £ {^,^}- Spin and isospin basis are taken in z ori- 
entation which means, for instance, that |\ff(x,t, 1T)| 2 is 
the probability of having at a; a quasiparticle with spin 
pointing along +z (f) of electron type (fr). 

Introducing a state wave number k we now assume the 
following expansion 



(fc) 



e lkx Xs 



{V<r)Xs T (Vr) , (3) 



where the quantum numbers are s CT = db, s T = ± and the 
spin/isospin two-component states fulfill 



TzXs t (Vt) = S T Xs T (Vr) ■ 

dk) 



(4) 
(5) 



The amplitudes ip s Js T in Eq. © can be viewed as a set of 
four complex numbers representing the state for a given 



wave number k. Substituting Eq. © in Eq. (fT]) we easily 
get the equations for the state complex amplitudes ip Str s T 
(we drop label k for convenience) 



h 2 k 2 

2 m 



-/i|s 

+A ^ S 



Abs ct + ak s a s T sm.(j>B — E 



Yk S a S T COS 4> B lAs7s T = , (6) 



where s = — s. Equation © is an algebraic (matrix) 
equation with nondiagonal couplings of ^ s<jSt to ips^ 
through pairing and to ips^s T through the Rashba inter- 
action. 

It is important to realize that Eq. (J6j) has solution only 
for restricted values of the wave number k. An allowed k 
which is purely real describes a propagating state while 
a complex k with nonvanishing imaginary part describes 
an evanescent state. Finding the physically allowed k's 
at a given energy E, even numerically, is a non trivial 
task. Notice, for instance, that the matrix in Eq. ^ is 
non Hermitian for an arbitrary complex k. While many 
numerical methods can deal with Hermitian matrices, the 
methods for non Hermitian operators are scarce. 

We have followed the idea suggested in Ref. [37j to find 
the allowed k's in a robust and numerically stable way. 
We arbitrarily choose a particular spin and isospin (s, t) 
and solve the following linear system 



[Eq.©] 



if (So 
if (S B 



= (M) 



(7) 



We have found that Eq. (J7J) is well behaved numerically 
and always has a solution for any wave number k. In or- 
der to find the physically acceptable k's we subsequently 
look for the zeros of the following function obtained from 
the solution of Eq. ©, 



l.h.s. of Eq. @ for (s a ,s T ) = (s,t) 



(8) 



Notice that the state amplitudes and wave numbers ful- 
filling Fio{k) = from Eqs. ([7]) and © are also solutions 
of Eq. © and, therefore, they are the physical solutions 
giving the complex band structure of the system. 

As an illustrative example, Fig. [T] shows a typical com- 
plex band structure. Notice that, contrary to common 
band representations, we choose the energy as horizon- 
tal axis since, within our approach, E is given and the 
wave numbers are inferred. We only display the positive 
energy solutions, the negative ones corresponding to the 
specular image with respect to E — 0. As in Ref. 0] we 
set our energy units using the Rashba interaction as a 
reference. Specifically, our length and energy units are 
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(9) 
(10) 



The propagating states in Fig. QJi agree with those al- 
ready discussed in Refs. @Hl|. The evanescent states 
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FIG. 1: (Color online) Band structure as a function of energy 
for A s = OAEso, A = 0.25£ SO , fJ, = O.lEso, 4>b = tt/4. 
Solid and dashed lines correspond to 5R(fc) and Sr(fc), respec- 
tively. Panel a) shows modes having only one nonvanishing 
component of the wave number: solid lines are for propagat- 
ing modes with Sy(fc) = while dashed lines are for evanescent 
modes with 5R(fe) = 0. Panel b) shows the complex fc's for 
evanescent modes with both real and imaginary parts. The 
units are defined in Eqs. ((9J and (|10[) . Thin gray lines are 
guides helping the eye to correlate the position of the band 
extrema. 



provide, however, a novel view. We distinguish two types 
of evanescent modes: those shown in Fig.[TJi have a purely 
imaginary wave number, while those with wave numbers 
having both real and imaginary parts are shown in panel 
b). Notice that propagating and evanescent bands can 
be smoothly connected by the extrema points. In panel 
a) this connection is clearly seen as horizontal parabolas 
of opposite curvatures, showing how a purely imaginary 
wave number transforms into a purely real one as the en- 
ergy increases. In panel b), of an evanescent state 
connects with the extremum of the corresponding prop- 
agating band, while the imaginary part disappears when 
entering the propagating band sector. For a given energy, 
there are always 8 modes, as corresponds to a quadratic 
kinetic energy and the 4x4 matrix in Eq. ©. 

Summarizing this section, we stress that the suggested 
approach is a robust numerical algorithm allowing the 
determination of the complex band structure for an ar- 
bitrary set of parameters. We will focus next on the 
description of edge states, in particular zero energy ones 
representing Majorana modes. 
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FIG. 2: (Color online) Densities of the Majorana mode (edge 
state at zero energy) corresponding to the complex band 
structure of Fig. [1] The definition of p SatST (x) is given in 
Eq. GU. 



B. Edge states 

The complex band structure contains the necessary in- 
formation on all possible eigenstates. In particular, lo- 
calized states for which the wave function fades away of 
some defect or interface are, in general, superpositions of 
evanescent states. Let us consider a semi-infinite ID sys- 
tem restricted to x > 0. A priori, we expect the abrupt 
edge at x — to allow the formation of localized states 
vanishing at x = and decaying towards x — > oo. Such 
state has to be a superposition of evanescent waves hav- 
ing 3(fc) > (where 5 stands for imaginary part). Let 
us assume that the set of such wave numbers is {fc}, con- 
taining N modes for a given energy E. The boundary 
condition an edge state has to fulfill is then 



e cue, 

ke{k} 







(11) 



.,0) 



where the CVs are complex numbers and the '4>s a 's T s are 
the state amplitudes defined in the preceding subsection. 
Equation (| 1 1 [) is a set of four linear equations, for 

(s a s T ) = (++, + — , — h, )j an( i N unknown CVs. We 

can immediately conclude that if JV < 4 no edge state 
is possible since there are less unknowns than equations 
in Eq. (TIT]) . That is, if there are less than 4 evanescent 
modes with SJ(fc) > only the trivial solution Ck = 
is possible, i.e., no physical solution is present. On the 
contrary, if N > 4 it will always be possible to find a set 
of nonvanishing C' k s. If N = 4 the nontrivial solution 
will be possible if the system matrix is singular (zero de- 
terminant). In that case the condition for the existence 
of an edge state is 



det{$*i r } = 



(12) 



For instance, in the case of Fig. Q] the requirement of 
at least 4 evanescent modes restricts the possibility of an 
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FIG. 3: (Color online) a) Complex band structure as a func- 
tion of A B for fixed E = 0, A = 0.25Kso, H = O.lSso, 
4>b = 0. As in Fig. [1] solid and dashed lines display real 
and imaginary parts of the wave number. The inset shows 4 
modes with a large value of 5ft(fc). b) Absolute value of the de- 
terminant Eq. (| 1 2 p . A vanishing determinant is the condition 
for the existence of an edge state. 



edge mode to energies E < O.lEso- We have also com- 
puted the determinant of Eq. (fT2|) finding that it vanishes 
only for E — 0, while it grows linearly for E > 0. There- 
fore, in this example a zero energy edge state is formed. 
This is a Majorana mode for which Fig. [5] displays the 
corresponding density profile, obtained as 



Pip) = Ps«s T {x) , 



.(*) 



(13) 



(14) 



As expected, the densities vanish at x — and decay 
with damped oscillations for increasing x. This behavior 
is in good agreement with previous results.— Notice also 
that because of symmetry the densities remain invariant 
after spin-isospin inversion. 

The above results can be taken as examples of the 
method and its capabilities and we now address the ques- 
tion of how to rationalize the existence of topological 
phases. It was advanced by Oreg et alZ that the sig- 
nature of an emergent topological phase is the closing 
and reopening of the k = gap in the band structure 
of propagating states as one Hamiltonian parameter is 
varied. Increasing the Zeeman energy in a parallel field 
geometry the critical parameter for the topological tran- 
sition is 
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FIG. 4: (Color online) Same as Fig. [3] for field orientation 
along y {<j> B = tt/2). 



such that zero-energy edge (Majorana) modes are present 

(c) 

only when Ab > A^ . In the transverse field geometry, 
no zero energy modes are ever present. Let us analyze 
this scenario from the evanescent band structure point of 
view. 

Figure [3J shows the evolution of the wave numbers at 
zero energy as a function of the Zeeman parameter. The 
field orientation is along x and, in this case, we always 
find 8 evanescent modes for any Ab- This means that 
there is no propagating mode for this energy. Notice that 
four evanescent modes with a rather structureless behav- 
ior and having a large real part of the wave number are 
shown in the inset. In panel b) we display the absolute 
value of the determinant Eq. (fT2")l as a function of A^. 
The determinant has a discontinuous behavior: it decays 
from its maximal value for As = 0; it has an accidental 
zero for A^ w O.lSEso, when all wave numbers in Fig. 
[5^, become purely imaginary; and it consistently vanishes 
for As > 0.27 Eso, after the point where ^s(k) vanishes 
for two of the modes. 

The results in Fig. [3] prove that Majorana modes exist 
above a critical field (0.27Eso m that case) which is in 
agreement with the A^ of Eq. ([TBI anticipated by Oreg 
and collaborators^ Not only this, the analysis in terms 
of evanescent modes also sheds light on the reason why 
the topological transition is signaled by the vanishing of 
a gap for propagating states. Indeed, when Ab = A^ in 
Fig. © the complete vanishing of k for two of the modes 
indicates that they loose the evanescent character and 
become propagating, but only for a single point in As 
and thus they are not seen in panel a) as fully developed 
bands. At this value of As, therefore, there is no energy 
gap for propagating states. Notice, finally, that Fig. [3J 
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also proves that the vanishing of the gap for propagating 
states occurs for a wave number k = 0. 

When the held is in transverse direction (y) , as in Fig. 
21 no Majorana modes are present. Quite remarkably, 
in this case there are propagating modes at zero energy 
(dark solid lines appearing for Ab > 0.25Eso)- The con- 
dition of at least four evanescent modes with ^s(k) > 
is fulfilled only when Ab < 0.25Eso but, even in this 
region, no edge state is allowed because of the non van- 
ishing determinant, as proved in panel b). 



III. THE 2D CASE 

We consider now the generalization of the complex 
band structure approach to include one extra spatial de- 
gree of freedom in the transverse direction to the wire. 
As before, we assume the wire to be infinite along x, but 
now motion along y is also possible with the restriction 
—Ly/2 < y < Ly/2. This corresponds to a hard wall con- 
finement in y within a square well of length L y . We call 
this a 2D case although in the literature the name quasi- 
1D is often used to emphasize that there is confinement 
in the transverse direction while along the longitudinal 
one motion is free. 



Hamiltonian 



The generalization of Hamiltonian ([T]) is 



\ 2m 2m 



AoT x + — (p x a y - Py<T x ) t z , 



\b& • n 



(16) 



where V(y) is the square well of length L y mentioned 



above. Besides y confinement Hamiltonian (1161) includes 
the transverse kinetic energy and a p^-dependent Rashba 
term that, following a usual convention, we call the 
Rashba mixing term. The mixing character of this term 
can be understood from the fact that the p y operator 
couples different square well eigenstates. 

The analysis proceeds in a similar way to the preceding 
ID case with the important difference that now the state 
amplitudes of Eq. become functions of y, 



(17) 



Accordingly, the algebraic Eq. © becomes a second order 
differential equation due to the dependence on p y and p y 



PI 



2 1,2 



K Z k 



2m 2m 



E 



±bs 



aksa-Sr sin 4>b 



-as T {ik s a cos 4> B - ^) fc r {y) = . (18) 



Another difference with respect to the ID case is the 
dependence on the magnetic field orientation. Now y 
and z transverse orientations are no longer equivalent. 
However, in Eq. (j!8[) we have again assumed a magnetic 
field in the xy plane with an azimuthal angle 4>b since, 
in this geometry, there are no orbital effects of the field. 

From the numerical point of view solving Eqs. (|18|l for 
an arbitrary k is much more demanding than the ID case. 
We have devised a practical algorithm^ introducing a 
uniform grid in y and an arbitrary matching point y m . 
The differential equation (|18|) on all grid points but y m 
is rewritten using finite difference formulas, with an im- 
portant caveat: crossing from left to right of the match- 
ing point is avoided using non centered finite difference 
formulas. Another condition is that the state amplitude 
vanishes at the edges of the y domain. At the matching 
point y m , instead of Eq. (|18[) . we require the conditions 



dy dy J Vs a s T Vym) 



1 , if (s CT ,s T ) = (s,t) , 
, if (s a ,s T ) yt (s,t), 

(19) 

where (s, t) are a pair of arbitrarily chosen spin and 
isospin labels. In Eq. (fK)|) we have used the notation 
dS L ' R > /dy to indicate derivatives at y m using only left 
(L) or right (R) neighboring grid points. Of course, a 
physically valid wave number must correspond to a con- 
tinuous first derivative of tp s t{y) at the matching point. 
We thus numerically determine k from the zeros of the 
function 



J°2D{k) - 



dy 



(20) 



As in the ID case, when J^dCO = is fulfilled the so- 
lution from Eq. (fH?)) is equivalent to that of the full Eq. 

The robustness of the suggested numerical algorithms, 
both for 2D and ID, is seen in the fact that for any com- 
plex k the T functions can be computed, avoiding singu- 
lar matrices or ill-behaved numerical problems^, Finding 
the physically allowed fc's is then accomplished scanning 
the complex k plane to determine the zeros of T with the 
desired accuracy. We have checked that the solution is 
not affected by the arbitrary choice of state components 
(s, t) and matching point y m . 



B. 2D edge states 

In 2D the boundary condition for an edge state with 
x > is similar to Eq. (TTTj) . with the state amplitudes 
changed from scalars to y-functions, 



E 

fce{fc} 



C k 4%(y)=0 



(21) 



We immediately notice that now the boundary condition 
requires an infinite number of evanescent states {k} due 



6 



to the infinite number of possible values of y. In practice, 
a truncation to a restricted set of N evanescent modes 
has to be done and, accordingly, the condition (|2"Tj) can 
be imposed on N/4 values of y (the number of modes 
should be a multiple of 4) . 

Instead of using specific values of y in Eq. (|21[) , we have 
projected Eq. (f2Tj) on the set of N evanescent modes, 
finding the matrix equation 

M k . k C k =Q, (22) 

ke{k] 

where 

M k , k =Y.Jdy 4 k > (y) 4% (y) . (23) 

When the number N of evanescent modes is large enough, 
the zero eigenvalues of matrix M. correspond to Majorana 
edge states. This Hermitian matrix can be diagonalized 
numerically with standard methods. For the ID case of 
the preceding section the diagonalization of matrix Ai is 
equivalent to the approach based on the determinant, Eq. 
(fl"2|) . since a zero eigenvalue is associated with a vanishing 
determinant. In 2D matrix M allows a truncation of the 
infinite number of evanescent modes. 

The complete determination of the evanescent bands 
in 2D is a very tedious task due to the large number of 
states when varying A^. Instead of finding the equiva- 
lent to Fig.|3]for the 2D case, we have followed a different 
approach to determine the existence of Majorana edge 
states. Since the topological phase transition is signaled 
by a complete vanishing of the complex wave number, we 
have scanned the value of Tiuik = 0) when varying the 
Zeeman parameter Ab of the Hamiltonian. The regions 
between nodes of T^d^ = 0) are the different phases, 
and we have just studied a few selected values of in 
each region. For a given A^, we look for a set of evanes- 
cent modes in the complex-/c plane. With these modes 
we then find matrix Ai k > k and compute its lower eigen- 
values. As discussed above, each vanishing eigenvalue is 
associated with a Majorana zero mode. 

In absence of Rashba mixing each transverse mode of 
the quantum wire behaves as an independent ID system: 
there is an onset for the appearance of the n-th band 
Majorana mode 

A<g] n = y/Al + (»-e n y, (24) 

where e n is the n-th eigenvalue of the transverse quan- 
tum well, and multiple Majorana modes are possible in- 
dependently. This physics is contained in the results 
shown as a dashed line in Fig. [5ji where, as could be 
expected, the derivative discontinuity J r 2c(fc = 0) has 
nodes precisely when A B = A^ n . We have also ob- 
tained that, in this case, matrix Ai has n different zero 
eigenvalues when evaluated for a value of Ab such that 
Ag n < Ab < Ag\ n+1 . In practice, a fast numerical 
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FIG. 5: (Color online) a) Value of the derivative discontinuity 
T2D for k = 0. Solid and dashed lines correspond to presence 
and absence of the Rashba mixing term, respectively. In ab- 
sence of Rashba mixing (dashes), the nodes signal the appear- 
ance of successive Majorana modes when increasing Ab ■ For 
the case with Rashba mixing, panels b) and c) show the evo- 
lution of the lower eigenvalues of matrix M when increasing 
the number of evanescent modes. The appearance of a zero 
eigenvalue is seen in b). Parameters: Ao = 0.25Eso, fJ> — 0, 
E = 0, L y = 5L SO - 



convergence of these lower eigenvalues to « 10 9 is eas- 
ily obtained. 

Conspicuous modifications to the above scenario are 
seen when including the Rashba mixing. Notice that in 
Fig. [5] we are assuming a regime of strong spin-orbit in 
which Eso ~ £ n +i~£ n - We find that for A B < 0.7lEso, 
below the first node of Fig. [SJi solid line, the matrix M. 
has no zero eigenvalue. In the region between the first 
and second nodes, Q.llEso < Ab < 2.VJEsoi we find 
clear evidence of a zero eigenvalue, as proved in Fig. [5Jd. 
Convergence is, however, rather slow and a large number 
of evanescent states is required. As a practical scheme, 
we have included the evanescent states in an ordered way, 
taking the modulus of the complex wave number as order- 
ing parameter. This criterion gives, for instance, that in 
Fig. Hb for I At I < 5.5i^o 

we include a total of 36 evanes- 
cent modes. 

In the region 2.17 'Eso < Ab < 3.97 Eso, between the 
second and third nodes of Fig. solid line, the lower 
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FIG. 6: Density distributions for the Majorana zero mode 
shown in Fig. [5b. Arbitrary units are used in a gray color 
scale. Contour lines are displayed for a better presentation. 

eigenvalues of matrix Ai converge to finite values, as 
shown in Fig. [5}; for As = 2.5Eso- This behavior is due 
to the Rashba mixing, for, in its absence, we find two zero 
eigenvalues in this region. Though numerical, the differ- 
ence between panels b) and c) of Fig.[5]is quite clear. We 
have also checked that for Ab = A.5Eso, above the third 
node, a similar behavior to that of Fig. [5t (without zero 
eigenvalues) is found. We conclude that the presence of 
the Rashba mixing hinders the coexistence of two or more 



Majorana modes, which is in agreement with the results 
of Ref. [30( for finite systems. It is remarkable that this 
hindrance persists even in the semi-infinite system. 

For completeness, Fig. |6] shows the density distribu- 
tions corresponding to the Majorana mode of Fig. [5b, 
using a 2D generalization of Eq. (JUJ). Notice that, on 
the scale of the figure, the density vanishes for x = and 
all values of y which is indicating the good convergence 
of the boundary condition with the number of evanescent 
modes. Of course, the densities also vanish asymptoti- 
cally for x — > oo and, as in the ID case, remain invariant 
when inverting spin and isospin. 

IV. CONCLUSIONS 

We have discussed the physics of Majorana modes in 
semi-infinite ID and 2D wires using the complex band 
structure approach. This formalism provides a natural 
description of the zero-energy edge modes as superposi- 
tions of multiple evanescent waves. The boundary con- 
dition of a vanishing wave function at the edge cannot 
always be fulfilled, this limitation defining the parameter 
regions (phases) where Majorana modes exist. 

The phase transition occurs when for zero energy the 
complex wave number of an evanescent band vanishes. 
The phase transition is seen as the emergence of a zero 
eigenvalue of matrix Ai, defined in Eq. (J23J) . In ID the 
analysis with matrix Ai is equivalent to the determinant 
of state amplitudes; while in 2D matrix M. allows to mon- 
itor the convergence of the lower eigenvalues when the 
number of evanescent modes is increased. In 2D we find 
evidence that, in case of a strong Rashba coupling, the 
Rashba mixing effect hinders the coexistence of multiple 
Majorana modes in the semi-infinite wire. 
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